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ABSTRACT 



We develop a formalism to investigate the behavior of quantum field and quantum 
Oh. ground state when the field is coupled to perturbation that periodically oscillates. 

Working in the Schrodinger picture of quantum field theory, we confirm that the 
phenomenon of parametric resonance in the classical theory implies an instability 
of quantum vacuum, and correspondingly it gives rise to catastrophic particle pro- 
duction if the oscillation lasts indefinitely; the produced number of particles expo- 
nentially increases without bound as time proceeds. The density matrix describing 
the limiting stage of the quantum state is determined by a small set of parameters. 
Moreover, the energy spectrum and the intensity of produced particles are worked 
out in greatest detail in the limit of weak coupling or small amplitude perturbation. 
In the case of strong coupling or large amplitude perturbation the leading adiabatic 
formula is derived. Application to cosmological fate of weakly interacting spinless 
fields (WISF) such as the invisible axion, the Polonyi, and the modular fields is 
discussed. Although very little effect is expected on the invisible axion, the Polonyi 
type field has a chance that it catastrophically decays at an early epoch without 
much production of entropy, provided that an intrinsic coupling is large enough. 



I Introduction 



What would be the fate of the quantum ground state if a periodic oscillation 
of some perturbation lasts for an indefinitely long time? This seemed a purely 
accademic problem, until modern unified theories of the micro-world predicted many 
weakly interacting spinless fields (WISF), and until it was clarified that condensate 
oscillation of these fields may become the major massive component of our universe, 
at least during some epoch of the early universe. Indeed, in some case (the invisible 
axion case U) this type of WISF becomes the cold dark matter candidate, even of the 
present universe, while in another case (the Polonyi field case [El) their decay creates 
too huge an entropy, destroying either the standard scenario of nucleo-synthesis or 
baryo-genesis, which then eliminates an interesting class of supergravity models. 

In this paper we shall formulate the fundamental problem of attacking the behav- 
ior of quantum field and quantum state under periodic perturbation. Our formalism 
is based on the Schrodinger picture of quantum field theory. In this formulation 
relationship between the quantum wave function and the classical oscillator motion 
becomes evidently clear. Thus we are able to confirm without any ambiguity that the 
exponential decay law of the initial ground state is related to the classical instability 
(parametric resonance) under periodic perturbation. Moreover, we are able to work 
out the density matrix of the quantum state initially started as the ground state, 
and able to compute the energy spectrum and the intensity of produced particles at 
any asymptotic future time. 

The phenomenon of catastrophic particle production that occurs at the limiting 
late stage of periodic perturbation has many features so distinct as to be verified 
unmistakably if it ever can be checked by laboratory experiment. It may also have 
important implication on cosmological WISF. We briefly examine the two cases of 
immediate interest, the invisible axion and the Polonyi field. It is shown that al- 
though very little effect is expected for the axion, the Polonyi type field has a chance 
of being dissipated away without much entropy creation. This would be a blessing to 
some phenomenology of supergravity models, but we leave to future the final word 
on this issue until more details are worked out. 

Some aspect of particle production under periodic perturbation has been inves- 
tigated in the past, especially in regard to the reheating problem after inflation 
@, ||. Aside from our different motivation on WISF, we point out that our main 
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achievement in the basic formalism is to clarify the nature of quantum state at the 
level of the density matrix, which has not been available so far even in the limited 
situation. Only with this level of understanding one can completely determine the 
limiting behavior of the quantum ground state, and can compute many quantities 
from the first principle. In clarifying the precise nature of quantum state, it is 
crucial to implement the unitarity, namely conservation of the probability. Our for- 
mulation manifestly satisfies the unitarity at every stage of approximation, whereas 
in some calculation in the literature the unitarity is obscured. We hope that we 
have contributed much to this regard. In technical aspects, we mention two major 
improvements here over previous works: explicit formulas from all higher resonance 
bands in the small amplitude or the weak coupling case, and derivation of a major 
adiabatic contribution in the large amplitude or the strong coupling case. 

The rest of this paper is organized as follows. In explaining the content of each 
section, we shall also state our main conclusion so that the reader can regard this 
as a summary of our work. In Section II we formulate the Schrodinger picture of 
the quantum ground state subjected to continuous periodic perturbation after an 
initial time. We only consider the simplest type of model of this sort, the Gaussian 
model with a periodic mass function. The wave functional of the entire system 
can be decomposed into a product of wave functions of Fourier modes since spatial 
translational invariance is kept intact. The Gaussian wave function of each mode 
is then written in terms of a single function of time u(t) that obeys the classical 
motion of field oscillator under the periodic potential. The width of the Gaussian 
peak is given by the modulus and the phase in the Gaussian exponent by the 

logarithmic derivative of the modulus. Thus analysis of the quantum wave function 
is reduced to that of the classical equation under some specific initial condition, 
determined from the initial condition of the quantum system. The spread of the 
Gaussian peak corresponds to excitation of higher levels of harmonic oscillator, thus 
production of particles. Hence growth of the field oscillator amplitude in the classical 
equation implies an increase of particle production in the quantum state. 

In Section III we explain how general Gaussian model is described by the density 
matrix in such a way that the unitarity is manifest. The density matrix in the 
Fock base is much easier, in the Gaussian model, to handle than the wave function 
itself, and is much more convenient to compute various observable quantities. The 
diagonal element of the density matrix, namely the probability distribution of n 
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particle production, is written in terms of a quantity that may be interpreted as the 
average number of produced particles. 

Section IV is a mathematical addendum in which we first explain how the in- 
stability region of the parametric resonance may be identified and derive the width 
of the resonance band in the case of small amplitude perturbation. Our method 
also makes it possible to compute the growing mode function in a systematic series 
expansion of the coupling factor. We then sketch how to obtain solution to the same 
classical equation under the large amplitude perturbation. The adiabatic approxi- 
mation is excellent in this regime. The classical amplitude of field oscillator is given 
in this case by an elliptic integral of the second kind, from which we determine the 
limiting behavior. This section contains some useful formulas and new derivation of 
old results that are not written in the standard textbook || on the Mathieu equation. 
Detailed derivation of weak coupling formulas is relegated to Appendix. 

In Section V we combine results of the previous sections in order to determine 
the limiting behavior of the evolved ground state. Under the small amplitude per- 
turbation there exist two vastly different time scales in the problem, the short time 
scale of oscillation, and the long time scale of decay or the inverse of the production 
rate of particles. We first time-average over the short time scale, ignoring transient 
phenomena. We then mode-sum over frequencies within the resonance band. In this 
way it becomes possible to obtain a useful formula of the spectrum intensity at any 
asymptotic time. Adding contribution from resonance frequencies in infinitely many 
bands gives rise to the exponential decay law of the ground state and its associated 
particle production. A precise relation between the two is established. If the periodic 
perturbation lasts indefinitely, the particle production becomes catastrophic. In the 
limit of small amplitude perturbation the spectrum of produced particles is equally 
spaced at to = |n (n = 1 , 2 , 3 • • •) of increasingly narrower band width, and the 
intensity of each band is determined as a function of time. In the large amplitude 
regime the particle production rate has a dependence on the oscillation amplitude 
quite different from the small amplitude case. This is shown by using the leading 
adiabatic formula in the previous section. 

Final Section VI describes cosmological application on WISF. We first explain 
briefly the cosmological damped oscillation of WISF. After the onset of the oscillation 
the Hubble time is time scale of the amplitude damping, hence during one Hubble 
time one may ignore effect of damping and assume undamped periodic oscillation. 
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After introducing a general parametrization on WISF properties, we discuss effect 
of the cosmological redshift that may push away relevant frequencies outside the 
resonance band. Rate of particle production while staying within the resonance 
band is a real measure of the decay rate of cosmologial WISF. It turns out that 
frequencies in general redshift away in time much shorter than one Hubble time 
under the small amplitude perturbation of WISF. This redshift effect is thus crucial 
in diminishing the particle production rate under the small amplitude perturbation. 
In most cases of cosmological application it is then important whether the large 
amplitude perturbation may or may not occur, since in this case the redshift effect 
plays a minor role. It is shown that the large amplitude perturbation is realized 
at least in one case examined: when a large intrinsic coupling is involved, which 
may or may not occur for the Polonyi field, depending on supergravity models. 
Using the formulas applicable under the large amplitude perturbation, we examine 
how the produced energy density relative to the parent's varies with cosmological 
evolution and under what condition this ratio becomes of order unity, in which case 
a catastrophic decay of WISF is expected. It is suggested that in the Polonyi field 
case the catastrophic particle production terminates the Polonyi field oscillation if 
an intrinsic coupling is larger than ~ 10. 



II Model in the Schrodinger picture 

In order to extract the most important features of the problem and not to be 
confused by minor details, we adopt a simple model of scalar boson field ip(x ,t) 
coupled to a classical oscillator £ (t) with the Hamiltonian, 

H = J d*x [ I + \ (Vy,) 2 + \ mV - \ l4(t)<f? ] • (1) 

A slightly general form of time dependence of the periodic oscillator £(£) is taken 
here: 

m = £+ e'^ + £_ e 2 * 2 , *=!*, (2) 

with real constants £±. The oscillation frequency is given by E and for a practical 
purpose one can take an amplitude of £+ = £_ = | , so that = £ cos(Et) . The 
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reason we take a more general form of £ + ^ is to show that none of our conclusion 
depends on the special cosine form, the result resting solely with presence of the two 
terms of the form, e ±2tz , in any finite mixture. On the other hand, a single term, 
either e 2lz or e~ 2lz , does not work. As mentioned in Introduction, we may envisage 
the situation in which the oscillator is a homogeneous component of a genuine 
field variable £(x,t). 

In the rest of discussion, we take the Schrodinger picture of the quantum Hilbert 
space. The field variable <p(x ,t) is decomposed into its spatial Fourier components 
qtif), and the wave functional of the state vector is written as 

*[{<K},*] = n^(ff S ,*)- (3) 

k 

Due to the space translational invariance it is obvious that each Fourier sector is 
linearly independent. For definiteness we enclose the system within a large normal- 
ization box of volume V such that each mode is labeled by a discrete index k. The 
wave function ijj^ obeys the Schrodinger equation, 

. di[> 1 d 2 ip 1 



+ (4) 



dt 2 dql 2 
with 

u 2 (t) = u 2 -^{t). 

To simplify our notation we denoted each Fourier component by the index uj = 
\j k 2 + m 2 , and we shall omit its suffix uj whenever no confusion may occur. 
The Schrodinger equation Eq.|] is solved by a Gaussian ansatz, 

^j{q u ,t) = -= exp( ~-ql), (5) 

y/U 2 U 

where it is found that 

d 2 u 



(]f2 + uJ 2 (t)u = 0. (6) 

The dot ' denotes the time derivative 4. This equation is equivalent to the corre- 
sponding classical field equation coupled to the oscillator £. Namely, the quantum 
problem is converted to the classical problem by the Gaussian ansatz. 
By changing the variables like 

u(z)=u(t), z= l -Et, h = (^) 2 , ± = ^, (7) 
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we find that the classical equation is of a generalized Mathieu type, 

d 2 u 



+ (h-9 + e- 2tz -9_ e 2tz )u = 0. 



dz 2 

Due to the reality of the constants 9±, once a special solution f(z) to the second 
order differential equation Eq.^ is found, another linearly independent solution is 
obtained as f*(—z). Furthermore, from general discussion || of the Mathieu type 
equation and as will also be explained in the next section, the special solution is 
written as the form of e Xz P{z) where P{z) is periodic: P(z + tt) = P(z) . The 
parameter A is either real or purely imaginary, indicating respectively unstable and 
bounded solution. In mathematics these are generally called Floquet solution, while 
in physics they are parametric resonance for A real and Bloch wave for A purely 
imaginary (with time replaced by a spatial coordinate). 

Our main interest lies in evolution of the initial ground state, 



^(g,0) = (^) 1 / 4 e-t-or ; Uo = . s j UJ 2_^ + + U ) (9) 

7T v 

hence we take as the initial condition to the classical differential equation, 

«(o) = (-r l/2 , -(o) = ^ . (io) 

71 U 

With this initial condition one may further achieve a simplification of the wave 
function. First, separating real and imaginary parts with u(z) = R(z) + il(z) , one 
finds that 

ii _ E R' + U' _ E R'R + I'I J'R-R'I 
u~l> R + iI ~ ~2 ^ R 2 + I 2 + 1 R 2 + I 2 J ' 

where the prime ' indicates the z— derivative 4- . By noting the constancy of the 
Wronskian of the two solutions, R{z) and I(z) , 

RI -RI = (—y 1 ^ = n, 
n 



„2 



one obtains 



ii 1 d , , l0 ^ii 7i 



SR- = _ _ In \u 



2 3- 



u 2 dt u \u\ 2 

1 q 2 i d 

thus, ip(q,t) = —f= exp[-— -{vr- ^-^W?}]- 



We have ignored a q— independent phase factor that do not affect physical observ- 
ables. This is a basic formula for our subsequent investigation. 
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The Gaussian wave function Eq.[TT| of the evolved ground state is characterized 
by the spread and a phase factor, the spread of the Gaussian peak being given by 
the modulus \u\ that almost, but not quite completely (rather u itself faithfully), 
follows the classical motion and the phase being given by the rate of time variation 
of this spread. Thus in our formalism there exists an unambiguous correspondance 
between the quantum wave function and the classical motion. Another nice feature 
of the formula EqjlT] is that the normalization of the wave function 



is automatically satisfied, ensuring the unitary development of the quantum state. 
Despite of the unitary evolution the ground state decays, and the decay amplitude 
is given by 



Thus exponential growth of the classical variable u in the instability band of the 
parametric resonance implies directly the exponential decay law of the initial quan- 
tum state. The decayed state is accompanied by production of field quanta, since the 
spread of the Gaussian peak implies excitation of higher harmonic oscillator levels. 
We shall have much more to say later on these issues after we introduce the density 
matrix of the quantum system. 



Ill Density matrix of Gaussian model and its general properties 

Although the wave function completely describes a quantum state, an equivalent 
description is possible in terms of the density matrix. It is sometimes, for instance 
in cases when observable quantities are needed, also useful to introduce the density 
matrix. In our subsequent discussion on the precise nature of the quantum state at 
late times, it is of crucial importance to use the density matrix formalism. This is 
due to that the density matrix in the Gaussian model has some simple properties, 
while preserving the fundamental relation of unitarity. 





(12) 
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The entire density matrix of our quantum system is decomposed into a product 
of submatrices p^ in Fourier modes as 

P = T[Pki Pk(Q ,<j',t)= 1>(q ,t)r(q',t) = (q\ > ( V> (*) W > ■ (13) 

fc 

Unless confusion occurs, we often omit the mode index uj and the time variable t. 
The wave function in our problem is of a Gaussian form, 

#z,i) = (— ^e"^, (14) 

7T 

with D(t) a known time-dependent function in our analysis of this section. All 
discussion that follows in this section is valid for an arbitrary D{t) provided that 
> 0. 

It is convenient to introduce the Fock base of the unperturbed field oscillators 

Qu>, 

n) = ^(alr\0), (q u \ n )=*(^g u +lj*)*ik(g u ), (15) 

Vn! V 2 V2cu oq^ 



ni 

2 



with ipo(qu,) oc e ^^l 2 the ground state wave function. The state vector in the Fock 
base is then expressed by 

(n\m) = (^) 1/4 ¥ J^ r dqH n (V^q)e-^^\ (16) 



n\ J-oo 

with H n (x) the Hermite polynomial of order n. The state vector components (n\ip(t)) 
trivially vanish for an odd integer n. In the even integer case, using an identity on 
the Hermite polynomial ||, 

r dye-y 2 H 2l (xy) = V^^r(x 2 -l) 1 , 

J— oo LI 

we obtain a very useful relation, 



dg Hn(^9) «" » ^ = ^ ^ t/"^ (^)' ■ (17) 
o l\ \ uo + D uj + D 

This formula greatly simplifies our subsequent discussion. 
The density matrix in the Fock base is thus given by 



(2/)! (2m)! ujMD |1/2 u - D t u - D* 

P-2I :2m 



l/2 U-^y U-^ 



2 l+m - 1 l\m\ \uj + D) 2 ' u + D' u + D* 
Pnm = (for norm odd) . (19) 

8 



It can be verified that this density matrix trivially meets the unitarity relation, 

oo 

Pnn = tr p = 1 , 

71=0 

as should be so by construction. The diagonal element of the sub-density matrix is 
written in terms of one parameter denoted by K, 



(21)! 2K .2-K « Auj^D 
P*>*-Z^\l2TK { 2TK } ' K = u* + \D\i- (20) 

A few immediate conclusions follow from consideration given here. The first and 
second moments of the density matrix give the average number and the dispersion 
of produced particles, which is computed as 

00 1 2 — K 

lN) = Yno nn = — —\uj-D\ 2 = - — — , (21) 

(AN 2 ) = (N 2 )-(N) 2 = 2\^^-\ 2 (N) 2 . (22) 

to — D 

The further D deviates from u, the more particles are produced. The diagonal 
element of the density matrix Eq.^0] is expressed solely in terms of the average 
number of produced particles ( N ) . 

The probability of the initial state being maintained is exactly given by 



\(0\m)\ 2 = Poo = ^^. (23) 

A precise relation of the decay law of the initially prepared state and the number of 
produced particles is then 

Poo= , = =■ (24) 

Thus when the initial state decays according to e~ r * , the produced number of 
particles increases like e 2Tt . This relation is exact irrespective of detailed form of 
the wave function, but it must be kept in mind that this exact relation holds only 
for a single frequency mode. We shall later derive a formula applicable when the 
mode-sum is performed. 
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IV Parametric resonance 



It is well known || that the Mathieu type equation with periodic coefficients 
exhibits instability in infinitely many band regions of the parameter space of h = 
(^) 2 and 6± = . In the weak coupling, or the small amplitude limit of h \9±\ 
the parametric resonance roughly occurs when the input frequency E is closely tuned, 
within a resonance band, to the characteristic frequency of the system, u, by the 
relation: lo = ~ n , namely , h = n 2 with n a positive integer. The parametric 
resonance occurs however small the coupling 8± is. Implication of this parametric 
resonance on the quantum ground state is the main theme of this paper. In this sec- 
tion we shall first review weak coupling analysis of the classical parametric resonance, 
explaining some details usually not written in the textbook on Mathieu's equation. 
We think that the material presented in the weak coupling limit either includes new 
results of great use in our later analysis, or is a new derivation of known results in a 
much more simple way than usually presented. In the second part of this section we 
sketch analysis in the large amplitude or the strong coupling limit. For brevity we 
only present results to the first leading term in the strong coupling limit. In the last 
part of this section we combine analysis in the above two regimes in order to identify 
strong coupling regions within a few lower resonance bands, which becomes useful 
in our later application. A precise identification of the adiabatic instability lines is 
difficult to find, by analytic means, in a general (6 , h) region, but for lower bands 
of n = 1 , 2 we can derive approximate analytic formula by comparing the stability 
chart obtained by numerical methods M to our weak coupling result. 

IV. A Small amplitude limit 

A Floquet type of solution is expanded as 

oo 

e Xz P(z)= c k (z) e l{n+2k)z , (25) 

k=— oo 

with the periodicity understood: P(z + tt) = P(z) . The basic idea of our approach 
is that within the narrow band near h = n 2 one may take e ±inz as the zero-th order 
approximate solution and assume small time variation of the amplitude Ck{z). It is 
thus legitimate to keep time derivative terms only in the two modes, k = and — n , 
in the form of co,-n ~ Aco,_ n , since |A| 1 . We anticipate a systematic hierarchy 
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of coefficients in the small parameter y/0 + 9_: 

Cfc = o[(M-) Min - (|fc|/2 ' |fc+n|/2) ], 

whose validity can consistently be checked by solutions derived. The growth rate 
factor A as a function of e = h — n 2 is obtained as described in Appendix: 



A < £ > = sV<ir> 2 -< £ - £ °> 2 ' (26) 



Here we defined various quantities by 



e o = o ~2 T e + e - ' for n ^ 1 ' 

2n -1 

= - - 0+0- , for n = 1 , (27) 

Aen ~ 2^-3[(n-l)!]2- (28) 
The case of n = must be treated separately. One finds that 



\ = J -^6+6--e , withe<-^M-- (29) 

The instability region with a real A is characterized by infinitely many narrow bands 
with the width Ae n Eq.^8[ However there is a gap between the the n = instability 
line Eq.|29| and the first n = 1 band: this region is defined by 

- ~M-<^<i-\/M^-~M-, (so) 

for a small 8 + 8_. In summary, the first few instability bands in the (9 , h) plane are, 
with 9 + = 9_=9, 

Q2 Q2 Q2 Q2 r 

h<--, l-9--<h<l + 9- J , 4--<h<4+-9 2 . (31) 

These formulas are accurate for small 9 values. At 9 as large as 3 the accuracy gets 
worse, but typically down to ~ 25%, which is still useful in many cases. 

In terms of the original oscillator variables the resonance band lies around uj = 
y n with the band width given by 

{^/e 2 t 

A ^ = 2^n[(n-l)!p ^ (32) 
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taking £ + = £_ = | £ for simplicity. The gap of stability region around the origin 
in the ( ^/JI^ , u; ) plane is specified by the three lines, 



cu = 0, ^ = 0, w = ^[l-2 



_ o Af€ _ l / 11/2 

2 [i m? 2 l m? jJ 



We call this region bounded by the three lines stability triangle. 

It is also possible to compute expansion coefficients of the solution u(z) in terms 
of a single overall factor c . We quote the leading result, 

qu/2 rjin/2 

u ~ c e Xz+inz + c„ n e Xz - inz = . = e **-H(n*-n) + + e \z-i(nz- v ) ^ 

^01 + 01 yf 9% + 02 

with 



tan 77 



1 



(n) a 

6+ " 6 4 n) = e o±^. (34) 



(n) 

e — e_ 



VI. B Large amplitude limit 

Having presented the necessary machinery in the weak coupling or the small 
amplitude limit, we proceed to the strong coupling or the large amplitude case. 
Here our discussion is limited to the leading approximation. Consider the Mathieu 
equation in a slightly different, but physically equivalent fashion: 

u" + Cj 2 (z)u = 0, Cu 2 (z) = h + 29cos(2z) , (35) 

with the initial condition, 

«(0) = *(„) =*>(.)). 

llX U 

We took both 9 > , /i > . Assumption of the initial ground state of harmonic 



oscillators means that cD(0) = v h + 26 1 is real and positive. For simplicity we take 
= such that 

The condition of the strong coupling or the large amplitude is 29 > h and we 
anticipate that this is realized by the adiabatic condition, 

(36) 



u 






> 1 u 1 


u 


LU 
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Neglecting the right-hand side of 

j- z [\WY + l -u*u>]=u>uu\ (37) 

that follows from the classical equation of motion, we find, by considering the initial 
condition, 

u(z) = 



, 1/2 exp[i f Z dzJ h + 26 - 40 sin 2 z ] . (38) 
2ix Jo 

The adiabatic condition Eq.^6| is obeyed for 9 ^> h as expected, except near the 
turning point of u(z) = 0. 

The exponent factor in the formula Eq.[J^ is an elliptic integral E(z ,k) of the 
second kind, 

iy/h + 29 [ Z dzJl-k 2 sin 2 z = Wh + 29 E(z , k) . (39) 
Jo 

with 

k 2 = , (40) 

29 + h ' K J 

which approaches 2 as 9 — > oo . The behavior of E(z , k) as a function of z for fc > 1 
is best studied by the complex integral in the cut w = sin z plane: 



/•sin z 

E(z ,k) — J dw 



1 — k 2 w 2 
1 — w 2 



Near the turning point w = ±-r one must analytically continue w into the complex 
upper half-plane so that one picks up the growing component which always dominates 
over the decaying component. Thus each time when z goes beyond the branch point 
singurality ± sin~ 1 i the exponent gets enhanced, and the total increment in passing 
through one period of Az = 2n is 



iAE(2n,k) = 4 dw\— — -. (41) 
Ji/k V 1 — w 

In the limit of 9 — > oo , namely k — > y/2 , this increment is numerically ~ 2.396. 
In the large time limit of z ^> 27r the growth factor accumulates, and one finds that 

/2^ e i.2V5?m i (42) 
V 9 E 
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When this formula is extended to an extreme, one may use it for the whole range 
of < uj < v7*? > although one must be careful of the stability triangle mentioned 
above. Furthermore, the important exponent factor has a variation at a finite hj (26). 
As a function of this variable the constant coefficient 1.2 in Eq.^ is modified to 

x/^ „/i 2 ■ + K f 1 j l2w 2 -l-5 . h 
J(<,) = 4(l--arcsm^— x _ w2 , 5=-. (43) 

The average value of this quantity for the range of 5 = — 1 is ~ 0.327, a factor 
~ 0.3 reduction from 1.2, when averaged with the phase space weight. We should 
consider this reduction for a better estimate. 

IV. C Strong coupling region within lower bands 

The large amplitude or the strong coupling formula above has been derived in 
the adiabatic approximation, with the condition h < 26 . But the actual region 
of applicable adiabatic region cannot be this simple, especially due to existence 
of infinitely many bands. From an inspection of the stability chart obtained by 
numerical methods 0] one can however gain a confidence of deriving approximate 
analytical expression of the instability region in which the adiabatic formula may be 
applied, at least for a few lower bands of n = 1 , 2 . 

We proceed as follows. First, we define the region bounded by the weak coupling 
band lines and the adiabatic line of h = 26; 

n 2 + e W < h< Min. (n 2 + e? } , 26) , e i n) ~ ± J^j ± ^ , (n + 1) . (44) 

Suppose then that 6 gradually increases from 0. For 6 < 0.33 the gap of the stability 
triangle exists and there is no adiabatic region. As 6 increases above 0.33, the 
adiabatic h = 26 line may cross the band boundary. Each time the adiabatic h = 26 
line crosses the lower weak coupling line h = n 2 + e_ , a new adiabatic region 
appears. By using the weak coupling formula Eq|JT] of the boundary line of the 
band, we obtain approximate expressions of the strong coupling regions within a few 
lower bands. These regions are classified according to the value of 6\ 

(1) 2 (V38 - 6) (~ 0.33) < 6 < 2 (V6 - 2) (~ 0.90) : 

q2 

1-6 <h<26, or 

8 

^[l- 2 4-^(4) 2 ] 1/2 <^< ™=lband(I), (45) 

2 mi 2 mi v 
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(2) 0.90 < 9 < 3 : 



n2 

< h <l + 6 - — , or 

8 

0<c<^[l + 24-i (4) 2 ] V2 : n = 1 band (II) » (46) 
2 2 

Q2 

for > 2 (4>/3 - 6) (~ 1.86) , 4 - — < /i< 26 , or 

^[4-^(4) ] V2 < <" < : n = 2 band (III) . (47) 



We call these three regions, Region I, Region II, Region III for later application. 
There is no simple analytic formula for 9 > 3 . But for a 9 larger than O[10] most 
of the region given by < h < 29 belongs to the instability region. 



V Limiting behavior of evolved ground state 

In this section we combine our previous formalism and mathematical analysis in 
order to determine the final stage of initially prepared ground state when the system 
is coupled to a periodic perturbation of frequency E. The sub-density matrix 
with uj = \j k 2 + m 2 of the quantum system is a bounded function of time for 
most of the modes, except within the infinitely many resonance bands, each joined 
to uj = -f n (n = 1 , 2 , 3 , • ■ ■) at the weak coupling limit of 9 + 9- = 0. Thus the 
limiting form of the density matrix as t — > oo is dominated by contribution from 
resonance frequencies. We therefore concentrate on combined effects from infinitely 
many bands, neglecting contribution outside the instability region. 

Recall, prior to any detailed discussion, that once the classical field oscillator 
amplitude u is given, non-vanishing sub-density matrix element p nm in the Fock 
base of frequency uj is obtained as 



_ v ( 2 ! (2m)! ujKD f uj-D uj-D* 

P2l,2m~ 2l+m _ ll[m[ Jim, Jim ~ \ ^ + ^ 2 \ ^ u + JJ > ^ u + JjJ > W 



7r i d , 1 

D =w + 2Jt x "W (49) 
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V.A Small amplitude case 



The leading small amplitude (9 <C h) behavior is computed from the previous 
general formula by puting in the initial condition. For simplicity we present results 
taking 6+ = 0- = ^§r ■ For the n— th band it reads: 

U (*)~i(^)-V2 
7T 

, 1 [ e inz {cos(2^) sinh(A^) - i sin(2??) cosh(A^)} + e~ inz sinh(A^) ] , (50) 
sin ( z/f) ) 

A = ^fE|En , e =(|)»-„>, ,„ M=2 /K^EI, (51) 

(^;/g 2 )" i g , a^„ . 

Ae n = — — , e± = ± 1 , 52 

2"-3[( n - 1)!]2 ' ± 2n 2 -l 2 vr; ' v; 

for e_ < e < e + with z — ^ t . Although it appears subleading, terms of order 
e -Az _ e ~\Et/2 mns i retained to get correct, non-divergent contribution near the 
threshold of the band, e ~ e± . 

As emphasized in the previous section, the density matrix element in each mode 
is governed by the two numbers, K ~ 1/(N) and a phase 5 with tan S = ^sD/u . 
Time dependence of this factor K~ 1 ~ (N) is necessarily complicated, because there 
are two different time scales involved, one scale very short of order 1/E and the other 
scale very long of order 1/(XE). Note that A = 0[ Ae ] < 1. A useful quantity is 
time average over the short time scale, for instance by averaging over a quasi-period 
of At = = ^ . Although exact time average of any density matrix element is in 
general difficult to compute in a closed form, it is relatively easy to get the average 
of the first moment, 

This was derived by using the technique of complex integration like 
i [ Z+ * d Z f(e^) = ± i — f(w) . 

TT Jz ^ ' 2lXl J\w\=l W ^ ' 

The number of produced particles thus computed has a correct finite behavior 
near the ends of the resonance band at e = e±, since the dangerous factor (e + — e) (e — 
e_) in sinh 2 (A2) and sin 2 (27/) cancel in the time average (N) . It is thus possible to 
get a finite result by keeping small 0[e~ XEt / 2 ] terms consistently. One indeed gets 
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a well-defined number of produced particles by summing over all unstable modes 
within the n— th band, 



u)~nE/2 



E 3 V Am 2 y/Ae^ exp[Ae n Et/(An)] 

■' n(n 2 - — — ) = , (54) 



2 8 v / 2vr 3 / 2 V ' E 2 ' ^Ei 

where we explicitly wrote n— dependence. An approximate estimate of integral, 

A_ sin z (2?7) Jo y J (At) 2 -Ay 2 4 V 2 

valid for Aez/(2n) 3> 1 was used here. We define the growth rate of produced 
particles by 

7ra = ^Ae n = \^! E2) \ m E. (55) 
r 4n 2™" 1 n[(n- 1)!] 2 v ; 

In the special case of n = 1 



TaTT" / 4m 2 / 9 „ , QE u£ 

Time dependence of the spectrum and the intensity of produced particles is of 
crucial importance in any application of the resonance production. The spectrum 
has an unmistakable feature of equal spacing, w = |?i, and of the band width, 

2(/</£ 2 r 

A ^-n[(n-l)!] 2 ^- 

Moreover, we explicitly gave the intensity rule at late times in Eq. [53]. If this phe- 
nomenon should be verifiable at laboratories, it has many experimental consequences 
that can be checked with detailed prediction. The total number of produced particles 
at any finite time, summed over all band index n, is finite, but becomes exponentially 
large at late times. 

The dominant contribution out of infinitely many bands is from the fundamental 
mode of n = 1, the case of pair production of equal energy -|. The dominance of the 
fundamental mode is true except in the large amplitude case of 9 = 2fi£/E 2 — 0[1] , 
for which many pairs of higher energy may readily be created. We shall have more 
to say on the large amplitude case in the second part of this section. Production 
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of high energy particles may superficially appear to violate the energy conservation, 



uj > 7^. The apparent violation is allowed due to that the external oscillation makes 
the concept of energy conservation meaningless at a formal level. However minor it 
may appear in the small amplitude case, it is important to keep in mind that due to 
their peculiar nature high energy particles may play important roles. 
The decay law of the initial vacuum is derived from 

IIPooH^nV^^^Pt-^E 111 ^)] « exp[-^A^t]= e - rw , (57) 

UJ UJ UJ ^ UJ 

which yields 



r = £r„, r. = W £ n/2 AM = i&S V" " -W 2* [(n - W ' (58) 

T is the total decay rate of the initial vacuum per unit volume. By taking the 
dominant n = 1 mode alone, one gets 

r -5.0 x l(r 3 (X) 2 , (59) 



for m <C E/2. The fact that the probability of the vacuum persistence Eq.|57] is 
suppressed by the volume factor V in the exponent is related to the homogeneity of 
space: the decay may occur anywhere with equal probability. The decay rate Eq.[55| 
may be compared with the Born rate of the usual one particle decay of mass M that 
conserves energy, 

{j 2 / Am 2 
1 ~ 32?rM V ~ ~M 2 ' 

This formula can also be derived from the above partial rate 



r t/ i 4m2 



by taking with E = M the amplitude of ^ = y jjy , which follows by equating 
the energy density of ^—oscillation | M 2 £ 2 to one particle mass density within the 
volume V, M/V . Difference here from this conventional Born formula lies in that 
the amplitude of oscillation £ may become much larger than that of one particle 
decay; indeed, the rate is enhanced by V ■ M£ 2 /2 which may become a huge factor 
in cosmological application. 
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V.B Large amplitude case 



We quote very briefly the main result. For simplicity we take £+ 




as before. First, the produced number of particles for each mode is given, using 
numerical values of the integral, by 



As explained in Section IV.B, the rate factor has a mode frequency dependence 
given by J(oj 2 /h£) with J{5) defined by Eq.|43|. The mode sum is performed with 
the instability constraint, which restricts the momentum range. If the whole phase 



adiabatic regime is restricted in more complicated ways as described in the previous 
section. Here we merely give an order of magnitude estimate by taking the full phase 
space and ignoring the u— dependence of the exponent. We postpone discussion of 
the actual reduction factor of the phase space to when application is made. Thus, 
we find for production of massless particles with m <C y/Jj£ , 




2.7(0) ~ 2.4. 



(60) 



space is available with the condition h < 26 alone, then k 2 



< /i£ — m 2 , but the 




(61) 



energy density 




(62) 



From 




the decay law of vacuum is 



Poo oc e 



-TVt 



with r 



(PO 2 , A, 



) ~ 2.4 x 10~ 2 (/z£) 2 , 



(63) 




again taking the full phase space. 
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VI Application: fate of cosmological WISF 



Modern unified theories predict many spinless bose fields that couple to ordinary 
matter extremely weakly. In cosmology based on these unified theories, there appears 
a grave potential problem. The condensate of these fields at the early universe may 
dominate the mass density of the universe after some epoch of oscillation, yet they 
may or may not decay due to their weak coupling. In either case they are likely to 
cause a disaster, because they either produce too much entropy when they decay 
at too late an epoch, destroying the standard nucleo-synthesis or the baryo-genesis 
scenario, or their mass density overdoses the present universe, which cannot be true. 
We shall examine how the new mechanism of resonant particle production discussed 
in previous sections affects this problem. In particular, we pay special attention to 
the case of the invisible axion [|lj, the Polonyi 0, and the modular fields, which may 
be of immediate interest. 

Weakly interacting spinless fields (WISF), which we generically denote by £, 
undergo damped oscillation at the early universe, obeying the evolution equation, 

£+3-£ + m^ = 0, (64) 

with a(t) the cosmological scale factor in the flat Robertson- Walker metric of ds 2 = 
dt 2 — a 2 {t) dx 2 . For m^t 1 , the amplitude of £ field behaves as 

m = £o [ ] 3/2 cos [ m « (* - *>) ] = £o ( £ ) 3/2 cos[m ? (t -to)}, (65) 
a[t) lo 

with to = 1/ rri£ , and To the temperature scale of the onset of damped oscillation. 
The initial temperature T is determined by equating the Hubble time scale to the 
inverse mass of £ field, 

^o = (^) 1/4 v^?^r, d=^N, (66) 

where N is an effective number of massless species contributing to the energy density 
at the temperature T . 

We consider a classs of WISF's that include the axion, the Polonyi, and the 
modular fields. Their properties relevant to our discussion are summarized by a 
few quantities, the mass m^, the initial amplitude of oscillation £o, and their generic 
Yukawa- type (| £y? 2 ) coupling denoted by /i, to a typical ordinary bose matter ip such 
as the photon, the Higgs, and the scaler fermion fields. To the Hubble time scale 
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the damping of oscillation amplitude is small, and one may assume exactly periodic 
oscillation. Mode functions of the bose field (p obey 

+ [u 2 -fi£(t)]q u = 0, u 2 = k 2 + m 2 , (67) 

during one Hubble time. Rewriting this equation in the standard form of Mathieu 
type gives the parameter set, 

E = mt, 0= 2 -^(yf 2 ^O o (^f 2 , h = A 2 = n 2 + e. (68) 

The initial 9 value, 9 = , is the most crucial quantity for the cosmological fate 
of f field. 

Effect of the cosmological redshift is very important in the case of small amplitude 
perturbation of 2\9\ < h, because it may readily push away resonance frequencies 
beyond the narrow band. If a frequency redshifts outside the band, the instability 
of this particular mode ceases to operate. Under the small amplitude perturbation 
the time At during which a particular frequency stays within the n— th band is 
determined by 

a(t) 



a(t + At) 
Assuming a small Ae n , one derives 



< Ae„ . (69) 



AtH < Ae n . (70) 

With Ae n <C 1 , the walk-out time is much smaller than the Hubble time if -1 . 

At the walk-out time Atn = ^ff- , the exponent factor in the produced number 
density is 

n, At {A6n)2 m * ^ ( T \^ 2 (7V\ 

Instil- ^ H ~ 2±n-A n ^ n _ l) ^y T J > W 

assuming the temperature dependence of the Hubble parameter in the radiation 
dominated epoch: H oc T 2 . Growth rate R n while staying within the band is then 
given by dividing the number density of produced particles by the number density of 
parent £ particles, \ £ 2 . By taking one species of a massless particle for simplicity, 
the production rate reads as 

r n3/2 r^ 2 ( T )~ 2 c^\r rV-'i r 16 ^' m W n (72) 
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The ratio of produced energy density to the parent's is 

^ = ^Rn. (73) 

Pt 2 

In particular, the dominant fundamental mode of n — 1 contributes with the rate, 
Rl * 4.0 x ID" 3 (^) 2 (^)- 2 exp[4 (4) 2 L ] . (74) 

This formula is valid only in the radiation dominated epoch. In the £— matter dom- 
inated epoch the temperature dependence in the exponent is modified to (T/T ) 3//2 . 
The prefactor (fr)~ 2 increases as time proceeds, hence as T decreases. This cannot 
be true for an indefinitely long time. Indeed, our asymptotic formula of produced 
number of particles is valid only for \ X n m^ t > 1 , that is, when the exponent above 
is large, hence when ^- is not too small. Here the strong coupling or the large 
amplitude problem is involved. 

The production rate under the large amplitude perturbation depends on various 
factors in different ways from the small amplitude case, in particular on the oscillation 
amplitude £, which may drastically change the fate of WISF. The frequency range 
that contributes to the growing particle production depends both on the initial 9 
value and on the temperature scale in question T, relative to T . For simplicity we 
shall only consider a moderate range of the initial parameter 6 = ^ir , for instance, 
O[20] or even less, which turns out to be sufficient for later application. The large 
amplitude perturbation region available is then restricted, initially to a wide h < 29 
region and at later stages to a few adiabatic regions, Region I, Region II, Region 
III, etc., as explained in Section IV. C. One must check in each case of application 
whether the necessary condition of applicability is fulfilled or not. Before going 
into any details of how this is met, we shall give an order of magnitude estimate of 
production rate in the adiabatic region. The number of produced particles of one 
light species of particles with m <C \fp£n (|^) 3 ^ 4 U P to one Hubble time is 



R = <^ ^ 2.0 x ID" 2 exp[ ( 2.4 - 0.65 ) 



0(^1,(75) 



^) 5/4 L (76) 



= 1.5 x 10- 2 (-^) 2 exp[ ( 2.4 - 0.65 ) . ... 

extending the phase space region to the full < uj < v / A^ • Two numbers written 
in the exponent correspond to the adiabatic approximation with and without the 
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to dependent reduction factor of the rate coefficient. The available phase space is 
actually reduced from this one, by a reduction factor of the form, f[^r (fj^) 3 ^ 2 ] • 
A precise form of the function f(x) will be given below when a specific application 
is made. The crucial exponent factor increases like (T0/T) 5 / 4 as time proceeds, in 
contrast to the T/T factor in the small amplitude case Eq.|74]. The precise temper- 
ature scale dependence (T /T) 5//4 here is trusted only when the coupling factor /x is 
independent of the particle energy uj, which does not hold in the axion case having 
a different (Tq/T) 1 / 2 behavior. 

The production rate given here is a yield produced during one Hubble time after 
the epoch of temperature T. There always exists accumulation effect of produced 
particles, but with a rapidly changing production rate the main contribution comes 
from one Hubble time at the maximal production. We therefore ignore a minor 
accumulation effect. 

With these in mind we proceed to discuss important cases of application. First, 
the invisible axion. The axion field which we denote by A has an effective two photon 
coupling in the form, 

C cS = —^AE-B. (77) 

J A 

The fine structure constant a appears, because this term is induced by a triangular 
loop of fermions. j a is the scale at which the Peccei-Quinn symmetry is broken and 
we ignored 0[1] factors in the coupling. The other coupling to matter fields such as 
quarks and leptons are irrelevant to our problem. The effective Lagrangian yields a 
modified Maxwell equation, 

•- - - a 1 • - 

E -V x B = --AB . (78) 

7T J A 

Relevant parameters for the resonant production is 

atrriA t , A 2 (fjn . 

V - - -r- w , 4o = J a , m A ~ — , (79) 

J A J A 

with A the QCD scale. The quantity #0 is dependent on the photon energy uj like 

e Q = 2 J^^^L. (go) 

m A 71 VTiA 



This factor is too small for ^- <ti it /(2a), and no appreciable resonant production oc- 
curs for the axion as far as the small amplitude formula Eq.[73] can be applied, because 
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frequencies immediately redshift away outside the band by the time Ae H 1 <C H 1 
I- 

Now the important question is whether there exists a region of parameters in 
which the large amplitude formula can be applied. This cannot occur in the axion 
case, because the stability gap constraint applied to the axion leads to incompatible 
conditions like 

-<--a(^) 3/2 , „>0.16m A [^(^) 3 / 2 ]- 1 . 

7T 1 7T 1 

Basically the presence of the stability triangle restricts 9 > 0.33 , which has no 
overlap with the adiabatic h < 26 condition due to the small - factor. Thus there is 
no way an appreciable resonant production takes place for the cosmological axion. 

The situation is different for fields of the Polonyi type. Here a generic dimen- 
sionless coupling may be of order unity or larger, but the effective coupling may be 
suppressed by a small ratio of an intermediate to a large mass scale. We express 
relevant parameters by 

F 2 F 4 
m?= M' ^ = M ' /i = ^°M^' (81) 

leading to the initial amplitude, 

#o = 22o- (82) 

The large mass scale M is the Planck scale and the intermediate mass scale F is 
of order 10 10 GeV or more in the case of the Polonyi field. We only consider the 
Yukawa coupling to spinless bosons, but there may exist coupling to the gauge field 
similar to the two photon decay of the axion, in which case the following estimate 
is modified in substantial ways. Assumption of the small amplitude oscillation leads 
to the growth rate while staying within the resonance band, 

^ = 0.99 x 10- 3 (|-) 4 (L)-* exp[4s 2 L ] . (83) 

Unless Qq is very large, for instance 250 or more, it is again unlikely that this type 
of WISF decays in the small amplitude regime. 

One must however use formulas relevant to the large amplitude perturbation if 
the exponent factor is large. The large amplitude formula yields production rate of 
particles in the mass range of m< v^o m C (f^) 3,74 > 

^ ~ 1.5 x 10- 2 gl (^) 4 /[ 2so(^r) 3/2 ] exp[ ( 2.4 - 0.65 ) y/gS (^) 5/4 ] , (84) 
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in the radiation dominated epoch. The two numbers in the exponent, 2.4 and 0.65, 
correspond to the adiabatic approximation with and without the u— dependent re- 
duction factor of the rate coefficient. A more refined estimate would lead to a 
modified exponent replacing the coefficient ( 2.4 — 0.65 ), which may weakly depend 
on g (T/T ) 3 / 2 . f(x) is the reduction factor of the phase space due to the restricted 
uo— range in the adiabatic region, which is given in regions of immediate interest, 

f(x) ~ f, (x>O[10]), 

1,5 5.-3 cc x , . 

= 1 -4<; +1 -M*><;- 1 + 2i>- ( 1J »<*< 3 ). 

= <k + \~^< (0.90<x<1.9), 

= <5 + £-£><§-S + I5>' (^<- 0.90). (85) 

Using these functional forms, one can readily find the temperature scale that the 
produced particle energy becomes comparable to the parent's. We did some sample 
computations. For instance, for g = 15 we find that the decay occurs at 

^~0.2, (86) 

taking F/M = 10~ 9 — 10~ 6 with the larger exponent 2.4, and 

T 

— ~ 0.06 - 0.09, (87) 
To 

with the smaller exponent 0.65. 

For a smaller g value one must take into account of the adiabatic constraint of 
6 > 0.33, since no sizable particle production is expected in the small amplitude 
regime. Ignoring the phase space suppression factor of order unity, the produced 
energy density roughly obeys 

-«l-5x 10- 2 9l exp [ ( 11 - 3.0 ) g^] , (88) 

Pi M 

right at 6 — 0.33 in the radiation dominated epoch. For this quantity to exceed 1 
so that the decay is completed in the large amplitude regime, the coupling g must 
satisfy g > 4 — 3 with the larger exponent 11, and g > 12 — 9 with the smaller 
exponent 3.0, taking for the intermediate mass scale F/M = 10~ 9 — 10~ 6 . Thus 
there seems no fundamental difficulty in getting a substantial decay of the Polonyi 
type field, provided that g > O[10]. The precise number of this bound needs to be 
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calculated by more sophisticated ways, but it should be clear that there is a definite 
possibility of resolving the Polonyi problem. 

One possible problem concerns effect of the Polonyi mass dominance on the par- 
ticle production rate. Modified time variation of the decayed energy density during 
the £— matter dominated epoch is given by 

P- = 1.5 x 10- gl (|-) 4 /[ 2, (|> 2 ] exp[ ( 2.4 - 0.65 ) ^ (^) 5/4 fe 3/4 ] ,(89) 

with T m the temperature at the onset of £— matter dominance, 

^ = |(— ) 2 - (90) 
T 3 m p i 

The rate in the matter dominated epoch increases more slowly than in the radiation 
dominated epoch, but still rapidly enough as the temperature T decreases. Right at 
the onset of matter dominance of T = T m the rate is given by 



M, 3 - 



m p i 



•exp[(2.4-0.65) v ^(A) 5 / 4 (^) 5/2 ]. (91) 
Taking M = m p i/\/8iT as an illustration gives the exponent factor as large as 

= iq ( 23-6.2 ) v /^ 

In order to complete the Polonyi decay at a temperature scale prior to the matter 
dominated epoch, there is a condition on £ parameters. These come from 9 > 0.33 
and ( ) > 1 , leading to 

9o ) 3 > 0.019, 92 

TOpi 

1.5 xl0- 2 s 2 (|-) 4 /[ 17£o (— ) 3 ] exp[(0.40-0.11) v ^( ! ^) 5/2 ]>l. (93) 
M m p i M 

Again taking M = m p i/^/8n as a guide and ignoring the phase space reduction 
factor of order unity, we find a constraint, 

1.5 x 10~ 2 gl {^f exp[ ( 23 - 6.2 ) ^] > 1 , (94) 

with g > 2.4. In the range of the intermediate mass scale of F/M = 10~ 9 — 10~ 6 , this 
constraint gives g > 13 — 6.1 for the larger exponent factor of 23 and go > 160 — 69 
for the smaller exponent factor of 6.2. 
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Although these numbers for the coupling seem a bit large, we note that the 
Polonyi field does not have to decay completely prior to the epoch of its dominance. 
Even in the £— matter dominated epoch the particle production rate increases, albeit 
with a slower rate than in the radiation dominated epoch. Right at the adiabatic 
limit of 9 = 0.33 the decayed energy density in the £— dominated epoch becomes 

-«l-5x 10- 2 gl (^) 4 exp [ ( 2.9 - 0.79 )g ^ ] , (95 ) 

again ignoring the phase space reduction factor of order unity. In order to get the 
Polonyi decay within the large amplitude regime, the right hand side quantity must 
exceed unity. With M = m p i/^/8ir as a guide, this happens for the larger exponent 
of 2.9 at g = 5.5 — 3.6, and for the smaller exponent of 0.79 at g = 19 — 13, taking 
the mass range of F/M = 10' 9 - 10~ 6 . 

A possible problem that might arise when the Polonyi field decays in the epoch 
of its dominance is a large entropy creation after the Polonyi energy is completely 
converted to radiation. This might dilute the baryon asymmetry created at the 
GUT epoch. We may estimate the entropy dilution factor S by assuming instanta- 
neous conversion of the Polonyi energy into thermal radiation. The reheated thermal 
temperature Tr is determined by 

d R T R = d m T^)\ (96) 

where the temperature T is obtainted by solving — — 1 . Ignoring a difference in the 
massless degrees of freedom di in different epoch, one finds for the entropy dilution 
factor, 

S~(^) 3 ~(^) 3/4 - (97) 

The further T deviates from T m , the more entropy is created, but with the dilution 
factor governed by a not too large exponent factor of 3/4. With M = m p \/ ^/8n the 
ratio of produced energy density to the parent's becomes 

j- = 1.5 x 10- 2 gl (^) 4 exp [ ( 23 - 6.2 ) ^ (^) 3/4 ] , (98) 

under the large amplitude constraint g (^-) 3 ^ 2 > 2.4 . One may thus compute the 
dilution factor S^go) for go > 2.4. For instance, with g = 10 — 50 the entropy 
dilution factor is computed as S ~ 4 — 2 at most, which may be regarded as a 



27 



minor effect. Given many uncertainties on the Polonyi field, especially on values 
of g and M , we believe the situation promissing for the resolution of the Polonyi 
problem. 

One great uncertainty not considered here is that the Polonyi mass may not be 
as simple as a constant m^. In general the Polonyi mass may depend on the field 
amplitude £. Furthermore, the Polonyi field may have decay channels other than 
spinless bosons, in particular massless gauge bosons like the photon and the gluon. 
In this case our estimate will be modified. For instance, if the Polonyi field decays 
dominantly into massless gauge boson pairs, then the exponential growth terminates 
at the onset of the Polonyi mass dominance in the case of go not large enough, 
and the ordinary decay at a much later epoch might create too much entropy. A 
more detailed analysis based on particular models is called for, but the catastrophic 
decay discussed here may potentially give a great impact on phenomenological model 
building of supergravity models. 
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Appendix Parametric resonance in the weak coupling 



We follow in spirit the physicist approach with guesswork and consistency check 
of results, by generalizing analysis of Ref.|§. A Floquet type of solution is thus 
expanded as 



e Xz P(z) 



£ Ck(z)e 



i(n+2k)z 



with the periodicity understood: P(z + 7r) = P(z) . Here we focus on the weak 
coupling or the small amplitude case, thus the instability region is expected to be in 
a narrow band. 

The basic idea of our approach is that within the narrow band near h = n 2 one 
may take e ±mz as the zero-th order approximate solution and assume small variation 
of the amplitude c&(z). It is thus legitimate to keep time derivative terms only in 
the two modes, k = and — n , in the form of c - n ~ Ac _ n , since |A| <C 1 . We 



anticipate a systematic hierarchy of coefficients in the small parameter a/# + #„: 



Ck 



>Min. ( |fc|/2,|fc+n|/2) ' 



With h = n 2 + e , |e| <C 1 , the original equation leads to an infinite set of relations 
among coefficients, 



7/c c k - 9+ c k+1 - 6L c fc _i = , 
7 fc = -4/c(n + fc) (k^0,-n), 70 



(99) 

e + 2mA, 7_ n = e - 2mA. (100) 



In matrix notation 

/0 




V 



7- 



\ 




oy 



/ • \ 

Ck 
C-k 

V • J 



(101) 



To proceed further, it is useful to separately treat the three blocks defined by 
k > 1 , k < — n — 1 , — n+l</c<— 1, and leave to the last the main Co (fc = 
0) , c_ n {k = —n) equation. Noting that 7^ = 7_ n _fc, one finds that the blocks 



29 



with k > 1 and k < — n — 1 yield identical results and 



D 



N 



V 



= 0-c 


lim (D N v )n 

N~*oo 


, c 


In 


-0_ 




-o + 


7at_i - 6L 




































72 








-0. 



(103) 



The central block of — n + l<A;<— 1 is then solved for c_ i and c_ n+i in terms 



of c and c_ n , yielding 

c_i = [( J B- 1 ) 11 + c o + ( J B- 1 ) 1 , n -i^ c_„ 



C-n+l = [(£' 1 )n-l,l 0+ Co + 1 ) 7l _i j „_i 0_ C. 



(104) 
(105) 



with i? = 



( l x -Q_ o 

-0+ 7 _ 2 -0_ 



o 









V 








(106) 



o -e + 7 _„ +2 -e_ 
o -0+ 7_„+i y 

Finally, we insert back c±\ and c_„±i into the last Co (/c = 0) , c_„ (A; = — n) 
equation. Noting that E^ 1 = E~\ n __ l , one finds that 

e + 2tn\-6 + 6_(D n 1 + E n 1 ) -OtE^ \( c 

- 0% E-_\ A e - 2mA - 0+0- (D u l + ^iV) ) { C- . 

Nontrivial solution exists only for 

(2nA) 2 = (0 + 0-) 2 E-_\ A E^ - [ e - + 0- (D^ + Etf) ] 2 = (e+ - e) (e - e_) , 

(108) 
(109) 



0. (107) 



e ± = + 0- (LTi 1 + ^n 1 ) ± M- ■ 

We then compute the inverse matrix elements to lowest non-trivial orders. We 
thus find that 



11 (Of) V 1 / 2 



(110) 
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Thus the band width of parametric resonance is given by 

)"/2 

22n-3(( n ri)!) 

The growth rate A takes a simple form, 



Ae = KU+U - J . (Ill) 



= 0, (113) 



The two special cases of n = and n = 1 must be dealt with separately. First, 
n = 1 case. By using 

e + 2mA-M-0D~ 1 )n - 0- W c ° 

-6+ 6-2in\-0+9-(D- 1 ) 11 ) \c-! 

we obtain 

e± = + 0_ % ± ^M- ^ - ^ 0+0- ± VM- , (114) 

with D = limjv_ >0 c Dn ■ 

Next, we deal with n = case. In this case there is no linear term of Co, and one 
must retain the higher order term of c . Moreover, there is no intermediate block, 
and the final coeffficient equation reads as 

-0+Ci + 7oCo-0_c_i = 0, 7o = A 2 + e, (115) 

leading to A = ]j - ^ 0+0_ - e , e < - J 0+0- ■ (H 6 ) 

Thus, there is no band of n — 0, and a gap extending to n = 1 band exists in the 
region, 

- ^0+0_ <h<l - y/0+0- - \ d +°- ■ ( 117 ) 

It is also possible to compute coefficients in terms of an overall factor c . We 
quote the result, 



— = C- r T /2 e-^\ tan^J^, (118) 

C- n 0+ V 6 ~ C - 

C- k = ( ^fci 1 0+ C + , Li 0- C n ) 

0f0g^ (0 + 0_) fe / 2 e-^ (fc!) 2 (0 + 0-) ( - fc)/2 ^ n< , <n _ n nip x 
n +0 n 2 2 ( fe ^ 1 ) [(k - 1)!] 2 2 2 («- fc -!) [(n- l)!] 2 J ' 1 " " J ' 1 J 

c ^(- 1 ) V M(n + t )! flt -^- '^^ 

c_»-fc ~ (- 1)* 22t H ( n + t)! g + c ~" ' < 120 ' 
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As promised, the expansion coefficients of terms e ± mz + tkz systematically contain 
suppression factors of #' fe ' , 6 = ^ \9 + 9 + \ . To leading order the functional form in 
the n— th band is then 

u ~ c e Xz+inz + c_ n e A2 - iri2 = ~ + + gAz-^-r,) _ ( m ) 
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